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I.  INTRODUCTION 


The  present  study  examines  the  computer  simulation  of  cw  HF 
chemical  laser  unstable  resonator  performance.  Previous  modeling  ,c  of 
the  Fabry-Perot  and  stable  resonator  performance  of  the  Hellos  CL  II, 
arc  driven,  subsonic,  continuous  wave  HF  chemical  laser  demonstrated 
that  the  Fabry-Perot  and  stable  resonator  models  needed  to  be  separately 
baselined  to  Fabry-Perot  and  stable  resonator  data,  respectively.  The 
separate  baselining  appeared  to  be  accounting  for  some  aspect  of  the 
resonator  modes  which  Is  not  Included  In  the  resonator  models.  Based  on 
the  results  from  the  Fabry-Perot  and  stable  resonator  modeling,  unstable 
resonator  modeling  was  examined  with  the  goal  of  determining  whether  a 
different  set  of  baseline  parameters  is  required  to  model  the  unstable 
resonator  data  even  though  the  wave  optics  unstable  resonator  model 
Includes  mode  structure. 

The  results  of  the  unstable  resonator  baselining  process  indicated 
that  one  set  of  baseline  parameters  that  allow  the  model  to  predict 
unstable  resonator  data  as  the  size  of  the  resonator  and  size  of  the 
laser  varied  could  not  be  found.  It  was,  however,  possible  to  obtain 
different  sets  of  baseline  parameters  for  Individual  resonator  sizes. 
The  best  fit  parameter  set  for  the  2  mm  slit  case  accurately  predicated 
the  total  power,  but  also  predicted  a  20  nsec  small  amplitude 
oscillation  that  was  not  seen  experimentally.  The  best  fit  parameter 
set  for  the  4  mm  slit  case  accurately  predicted  the  total  power  and 
period  of  the  .time-dependent  oscillations,  but  underpredicted  the 
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amplitude  modulation.  The  best  fit  parameter  set  for  the  5  mm  slit  case 
was  a  set  of  parameters  obtained  from  baselining  the  Fabry-Perot  model 
to  stable  resonator  data^‘4.  This  parameter  set,  used  in  previous 
unstable  resonator  modeling,  accurately  predicted  the  period  and 
amplitude  of  the  time-dependent  oscillations,  but  overpredicted  the 
total  power.  Neither  the  2  mm  slit  nor  the  4  mm  slit  baseline 
parameters  predicted  scale  effects  as  the  size  of  the  laser  was  reduced 
by  a  factor  of  two. 

The  inability  to  obtain  one  set  of  baseline  parameters  for  the 
unstable  resonator  model  demonstrated  that  the  model  is  inadequate.  The 
limitations  arise  mainly  from  the  treatment  of  the  fluid  dynamic 
mixing.  The  limitations  imposed  by  the  mixing  model  did  not  prevent  the 
Fabry-Perot  and  stable  resonator  models  from  being  baselined  to  the 
data.  This  suggested  two  possible  explanations.  Since  the  actual  fluid 
dynamics  is  independent  of  the  resonator  and  the  unstable  resonator 
model  is  a  more  exact  treatment  of  the  optics,  the  wave  optics  model  of 
the  unstable  resonator  may  have  brought  out  inadequacies  in  the  mixing 
model  that  were  not  apparent  with  geometric  optics  models.  Another 
possible  explanation  is  that  the  way  in  which  the  optical  fields  sample 
the  media  in  an  unstable  resonator  has  brought  out  the  inadequacies  in 
the  mixing  model.  Fabry-Perot  and  stable  resonators  have  no  upstream/ 
downstream  coupling  and  full  upstream/downstream  coupling  of  the  media 
and  optical  fields,  respectively.  The  unstable  resonator  couples 
upstream  of  the  optical  axis  and  downstream  of  the  optical  axis,  but  not 
across  the  optical  axis.  Both  possible  explanations  demonstrate  that 
unstable  resonator  data  as  a  function  of  resonator  size  and  size  of  the 
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laser  provide  a  severe  test  of  a  coupled,  wave  optics,  chemical  kinetic, 
fluid  dynamic  model  of  a  chemical  laser. 

The  unstable  resonator  model  was  used  with  the  best  fit  parameter 
sets  to  study  the  time-dependent  oscillations  in  an  unstable 
resonator.  Computer  calculations  indicated  a  correlation  between  the 
number  of  passes  for  a  wave  to  exit  the  resonator  after  leaving  the 
central  Fresnel  zone  and  the  oscillation  period.  This  correlation 
suggested  that  the  time-dependent  oscillations  are  the  result  of  a 
competition  between  chemical  pumping  and  radiative  deactivation  of  the 
upper  laser  levels;  the  oscillations  occur  only  if  the  media  is  not 
strongly  coupled  to  the  optical  fields  diffractively  or  geometrically. 
Computer  calculations  with  the  unstable  resonator  model  supported  the 
proposed  mechanism.  The  calculations  showed  that  the  oscillations  were 
damped  by  Increasing  the  rate  of  the  chemical  pumping  reactions,  and  the 
period  was  increased  by  increasing  the  degree  of  coupling  between  the 
media  and  the  optical  fields. 

In  Section  II,  unstable  resonator  computer  modeling  is  examined  and 
baseline  parameters  are  determined  for  different  resonator  sizes.  Ir 
Section  III,  a  mechanism  for  the  time-dependent  oscillations  is 
proposed.  In  Section  IV,  several  concluding  remarks  and  suggestions  for 
future  studies  are  presented. 
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II.  COMPUTER  SIMULATION  OF  CHEMICAL  LASER 
UNSTABLE  RESONATOR  PERFORMANCE 

Computer  simulation  of  chemical  laser  performance  must  adequately 
model  three  coupled  elements  of  chemical  laser  physics.  These  elements 
are  the  fluid  dynamics  of  the  flow  field,  including  the  fluid  dynamic 
mixing  of  the  fuel  and  oxidizer  streams,  the  chemical  kinetics  occurring 
In  the  fluid  media,  and  the  optical  resonator.  The  computer  simulation 
of  the  Helios  CL  II  and  CL  I,  arc  driven,  subsonic,  continuous  wave  HF 
chemical  lasers  is  accomplished  with  the  computer  codes  BLAZE  II8,8, 

MN0R037,  MN0R03SR1,8  and  MN0R03UR8.  The  MN0R03,  MN0R03SR,  and  MN0R03UR 
codes  are  rotational  nonequilibrium  models  with  Fabry-Perot,  stable 
resonator  geometric  optics,  and  unstable  resonator  wave  optics  models, 
respectively.  The  BLAZE  II  code  is  a  rotational  equilibrium,  Fabry- 
Perot  model  which  calculates  all  fluid  dynamic  properties  from  input 
Initial  conditions.  The  BLAZE  II  code  is  well  suited  for  the  simulation 
of  laser  performance  in  terms  of  total  outcoupled  power,  power  contained 
in  each  vibrational  band,  and  the  length  of  the  lasing  zone.  However, 
since  8LAZE  II  employs  rotational  equilibrium  kinetics,  it  cannot 
compute  realistic  power  spectral  distributions. 

Because  BLAZE  II  is  a  very  inexpensive  code  to  run  when  simulating 
chemical  laser  performance,  BLAZE  II  calculations  are  always  performed 
first.  After  suitable  agreement  between  BLAZE  II  calculations  and  the 
data  Is  obtained,  the  rotational  nonequilibrium  codes  (MN0R03,  MN0R03SR, 

MN0R03UR)  are  used.  In  the  rotational  nonequilibrium  codes,  profiles  s 
for  the  velocity,  temperature,  and  pressure  of  the  flow,  the  mass  flow  " 
remaining  In  the  primary  stream,  the  mass  flow  remaining  in  the 
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secondary  stream,  and  the  ratio  of  the  thickness  of  the  mixed  stream  to 
the  geometrical  thickness  of  the  flow  are  input  data.  These  profiles 
are  obtained  from  a  rotational  equilibrium  model  such  as  BLAZE  II. 

The  flow  field  In  the  Hellos  CL  II  and  CL  I  laser  bodies  is  a 
complex  three-dimensional  flow^,  Fig.  1.  A  sketch  of  the  appearance  of 
the  H£  jets  is  shown  on  the  left  side  of  Fig.  1.  These  jets  are  well 
defined  with  the  chemical  reaction  occurring  at  the  interface  of  the 
primary  (oxidizer)  stream  and  the  H2  jets.  The  appearance  of  these  jets 
as  seen  looking  upstream  is'  shown  as  Section  AA  in  Fig.  1.  To  treat 
this  complicated  three-dimensional  mixing,  the  flow  was  idealized  as  a 
series  of  alternating  primary  and  secondary  (fuel)  slit  nozzles  as 
shown  on  the  right  side  of  Fig.  1.  The  mixing  calculation  with  BLAZE  II 
was  performed  utilizing  a  linear  scheduled  mixing  option^’®. 

The  BLAZE  II  code  integrates  the  governing  ordinary  differential 
equations  describing  the  chemical  kinetics  and  fluid  dynamics.  The 
integration  is  performed  in  the  flow  direction,  starting  at  the  nozzle 
exit  plane.  For  the  Helios  CL  II  and  CL  I  lasers,  the  nozzle  exit  plane 
coincides  with  the  point  of  H2  injection.  The  integration  requires  the 
specification  of  all  of  the  dependent  variables  in  the  equation  set  at 
the  nozzle  exit  plane.  For  the  fluid  dynamic  equations,  the  dependent 
variables  that  must  be  specified  consist  of  the  static  temperature  and 
pressure  at  the  nozzle  exit  plane,  the  cross  sectional  area  of  the  flow 
channel,  the  mass  flow  rates  of  the  primary  and  secondary  streams,  and 
the  chemical  composition  of  the  primary  and  secondary  streams.  To 
obtain  values  for  these  starting  conditions,  the  static  pressure  is 
measured  5  mm  downstream  of  the  H2  Injectors  and  the  static  pressure  at 


the  nozzle  exit  plane  (H2  injectors)  is  varied  until  the  pressure  at 
5  mm  downstream  agrees  with  the  measured  value.  The  cross  sectional 
area  of  the  flow  channel  and  the  mass  flow  rates  are  measured.  The 
static  temperature  of  the  flow  at  the  nozzle  exit  plane  (H2  injectors) 
and  the  mass  fractions  of  F,  O2,  SFg,  SFX  (x  =  1-5),  and  S02  are 

unknown.  The  mass  fractions  are  unknown  because  the  degree  to  which  SFg 

is  dissociated  in  the  arc  is  unknown.  In  addition  to  these  quantities, 
the  lengths  after  which  the  primary  and  secondary  streams  are  completely 
mixed  must  be  specified  as  input  data  since  a  scheduled  mixing  model  is 
used.  Thus,  to  model  the  Helios  CL  II  and  CL  I  lasers,  there  are  four 
a  priori  undetermined  quantities  that  must  be  specified.  The  four 
quantities  are  the  percent  SFg  dissociation  in  the  arc,  the  static 
temperature  of  the  flow  at  the  H2  injectors,  the  primary  mixing  length, 
and  the  secondary  mixing  length.  A  set  of  these  four  quantities,  called 
baseline  parameters,  must  be  chosen  so  that  the  model  will  accurately 
predict  laser  performance.  The  process  of  determining  the  values  of 

these  parameters  is  called  baselining. 

An  objective  of  this  study  was  to  examine  unstable  resonator 

modeling  in  light  of  Fabry-Perot  and  stable  resonator  modeling  results 

which  showed  that  both  the  Fabry-Perot  and  stable  resonator  models 

needed  to  be  separately  baselined  to  Fabry-Perot  and  stable  resonator 

1  ? 

data,  respectively  *  .  The  Fabry-Perot  model  was  baselined  to  Fabry- 
Perot  data  because  the  model  did  net  accurately  predict  the  data  when  a 
set  of  parameters  obtained  from  baselining  the  Fabry-Perot  model  to 
stable  resonator  data  was  used.  When  the  baseline  parameters  determined 
from  baselining  Fabry-Perot  model  to  Fabry-Perot  data  were  used  in  the 


stable  resonator  model,  the  model  did  not  accurately  predict  stable 

resonator  data.  Thus  the  stable  resonator  model  was  baselined  to  stable 

resonator  data.  The  separate  baselining  of  the  Fabry-Perot  and  stable 

resonator  models  appears  to  be  accounting  for  some  aspect  of  the 
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resonator  modes  which  is  not  included  in  the  resonator  models  *  .  The 
goal  in  this  study  was  to  determine  whether  a  different  set  of  the  four 
baseline  parameters  is  required  to  model  the  unstable  resonator  data 
even  though  the  wave  optics  unstable  resonator  model  includes  mode 
structure.  To  serve  as  an  aid  in  discussing  the  unstable  resonator 
modeling,  brief  summaries  of  the  unstable  resonator  model,  previous 
unstable  resonator  calculations,  and  the  unstable  resonator  data  are 
given. 

2.1  THE  UNSTABLE  RESONATOR  MODEL 

The  efficient,  rotational  nonequilibrium,  chemical  kinetic,  fluid 
dynamic  model  denoted  MN0R037  was  coupled  to  the  Air  Force  Weapons 

Q 

Laboratory  physical  optics  strip  resonator  model  .  The  resulting  code 
is  denoted  MN0R03UR.  The  geometry  of  the  confocal  unstable  resonator  is 
shown  in  Fig.  2.  The  MN0R03  model  assumes  that  the  flow  and  chemistry 
are  independent  of  the  transverse  coordinate  z.  The  optical  cavity  is 
modeled  as  a  strip,  confocal  unstable  resonator.  The  beam  is 
diffractively  propagated  through  the  cavity  using  a  solution  of  the 
Huygens-Fresnel  integral  which  is  obtained  by  a  fast  Fourier  transform 
technique.  The  optics  model  is  coupled  to  the  chemical  kinetics  model 
through  the  thin  skin  gain  approximation  in  which  the  active  medium  is 
assumed  to  be  a  thin  sheet  in  front  of  the  large  mirror.  The  solution 


FLOW  DIRECTION 


MIRROR 


Figure  2  Confocal  unstable  resonator  geometry. 
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is  an  iterative  one.  The  calculation  starts  by  passing  a  plane  wave  of 
about  the  expected  intensity  (obtained  from  a  rotational  nonequilibrium 
Fabry-Perot  model  such  as  MN0R03)  through  the  empty  cavity.  The 
resulting  intensity  distribution  is  sent  to  the  chemistry.  The  species 
equations  are  then  solved  with  the  given  intensity  distribution  and  the 
thickness  of  the  mixed  flow,  Le,  and  the  gains  on  each  spectral  line  are 
calculated  and  sent  to  the  optics.  The  stored  wave  from  the  preceding 
optics  calculation  is  then  propagated  through  the  cavity  and  incremented 
by  the  gain  distribution  from  the  chemistry.  This  process  is  repeated 
until  the  difference  between  successive  iterates  is  less  than  some 
prescribed  amount.  Because  each  iterate  in  the  solution  of  the  steady 
state  equations  corresponds  to  one  round  trip  of  the  radiation  through 
the  cavity,  each  Iterate  may  be  regarded  as  a  time  step  in  the 
development  of  the  steady  state  solution10,11,12.  Thus,  fluctuations  in 
the  power  from  iterate  to  iterate  may  be  regarded  as  time-dependent 
oscillations. 

The  period  of  the  time-dependent  oscillations  is  calculated  by 
multiplying  the  number  of  iterates  in  one  period  by  the  photon  round 
trip  transit  time.  For  a  mirror  spacing  of  100  cm,  the  round  trip 
transit  time  for  a  photon  is  6.67  nsec.  The  oscillation  period,  then. 
Is  the  number  of  Iterates  for  one  period  multiplied  by  6.67  nsec.  The 
percent  amplitude  modulation  is  calculated  by  dividing  the  peak-to- 
valley  power  difference  by  the  average  power. 

All  of  the  Hellos  CL  II  and  CL  I  unstable  resonator  data  to  which 
the  MN0R03UR  model  was  compared  were  taken  with  Brewster  windows 
attached  to  the  laser.  There  were  no  vacuum  mirror  mount  data.  To 


Improve  the  comparisons  between  the  external  mirror  mount  data  and  the 
predictions  of  MN0R03UR,  Brewster  window  losses  were  added  to  the 
unstable  resonator  model. 

Cavity  losses  Introduced  by  the  Brewster  windows  were  modeled  in  a 
previous  stable  resonator  study*.  The  Brewster  window  losses  were 
determined  for  the  stable  resonator  model  after  the  model  had  been 
baselined  to  vacuum  mirror  mount  data.  In  the  stable  resonator  model, 
each  time  the  intensity  distribution  was  passed  through  a  Brewster 
window,  it  was  multiplied  by  an  effective  transmissivity, 

1  after  B.w!x^  =  *w  ^ncident^  (2-1"1) 

where  I(x)  Is  the  intensity  distribution  and  tw  is  the  effective 
transmissivity  of  the  Brewster  window.  By  varying  the  value  of  tw  and 
comparing  the  calculated  total  powers  to  stable  resonator  external 
mirror  mount  data,  a  value  of  the  effective  transmissivity  of  the 
Brewster  windows  was  determined.  For  both  the  left  and  right  Brewster 
windows,  tw  was  set  at  0.97.* 

Unlike  the  stable  resonator  model,  the  unstable  resonator  model  is 
a  wave  optics  model.  Thus,  the  unstable  resonator  model  performs 
calculations  in  terms  of  the  electric  field,  with  intensities  calculated 
from  the  electric  field  as  needed.  Since  intensity  is  proportional  to 
the  square  of  the  electric  field  strength. 


the  incident  electric  field  must  be  multiplied  by  the  square  root  of  the 
effective  transmissivity  used  in  the  stable  resonator  model, 

£  ..  _  , =  /F  £,  .  (2.1-3) 

after  B.W.  w  incident 

With  this  modification,  Brewster  window  losses  were  included  in  the 
unstable  resonator  model.  The  placement  of  the  Brewster  windows  in  the 
optical  cavity  is  shown  in  Fig.  2. 

The  value  of  tw  used  in  the  unstable  resonator  model  could  not  be 
determined  as  in  the  stable  resonator  case  since  there  were  no  unstable 
resonator  vacuum  mirror  mount  data  available.  Instead,  unstable 
resonator  calculations  were  performed  to  evaluate  the  use  of  the  stable 
resonator  value  of  tw  In  the  unstable  resonator  model.  The  calculations 
were  performed  with  tw  *  1.0  (no  losses)  and  tw  =  0.97  for  both  the  left 
and  right  Brewster  windows  for  the  low  pressure  Run  34  flow  rates.  The 
calculations  used  a  set  of  baseline  parameters  determined  by  baselining 
the  Fabry-Perot  model,  MN0R03,  to  stable  resonator  data^’4.  These 
parameters,  denoted  the  FPSR  parameters,  were  used  because  they  had  been 
employed  In  the  past  unstable  resonator  calculations®’^’^.  The 
unstable  resonator  calculations  are  compared  to  corresponding  stable 
resonator  calculations  in  Table  1.  Since  the  value  of  tw  from  the 
stable  resonator  model  resulted  in  total  power  reductions  that  were 
comparable  for  both  the  unstable  and  stable  resonator  models,  the  value 
of  ^  =  0.97  for  both  the  left  and  right  Brewster  windows  determined  for 
the  stable  resonator  model  could  be  used  with  the  unstable  resonator 


Resonator 


Total  Power  a% 

(Watts) _ 


Stable  resonator  data,  EMM* 


28.2 


MN0R03SR1  (tw  =  1.0) 

i 

|  MN0R03SR1  (tw  =  0.97) 

I 

Unstable  resonator  5  mm  slit  data,  EMM*^»** 


35.23 

25.55 

2.25 


27% 


MN0R03UR,  5  mm  slit  (tw  =  1.0)10’11 
MN0R03UR,  5  mm  slit  (tw  =  0.97) 


11.00 

7.11 


35X 


w 

V 

I. 

Table  1  Comparison  of  the  effect  of  Brewster  window  losses  on  stable 

and  unstable  resonator  calculations  for  the  low  pressure  Run  34 
flow  rates.  The  stable  resonator  baseline  parameters  were  used 
In  the  MN0R03SR  calculations  and  the  FPSR  parameters  were  used 
In  the  MN0R03UR  calculations. 
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2.2  PREVIOUS  UNSTABLE  RESONATOR  CALCULATIONS 

An  a  priori  calculation  with  the  unstable  resonator  model, 
HN0R03UR,  predicted  the  existence  of  time-dependent  oscillations  on 
lines  whose  saturated  gain  did  not  fill  the  unstable  resonator®' 

The  results  of  this  calculation  for  a  50%  geometric  outcoupled,  confocal 
unstable  resonator  with  a  large  mirror  diameter  of  1.0  cm  are  presented 
in  Figs.  3  through  5  and  Table  2.  The  calculation  was  performed  for  the 
low  SFg  and  low  H2  flow  rates.  Run  34,  at  a  laser  cavity  pressure  of 
5.4  Torr.  The  FPSR  baseline  parameters  were  used  in  the  calculation. 
There  were  no  Brewster  window  losses  included  and  thus  the  total  power 
was  overpredicted.  As  shown  In  Figs.  3  through  5  and  Table  2,  the 
calculation  predicted  a  40  nsec  oscillation  of  the  total  power  with  an 
amplitude  modulation  of  36%.  These  oscillations  occurred  on  lines  whose 
saturated  gain  did  not  fill  the  resonator.  On  the  Individual  lines,  the 
period  of  the  oscillations  was  approximately  40  nsec  and  the  amplitude 
modulation  varied  from  18%  to  100%.  Since  these  oscillations  were  not  a 
numerical  instabi 1 itylc  and  are  undesirable  in  a  laser  device,  a 
laboratory  experiment  was  devised  to  determine  whether  the  oscillations 
occurred  as  predicted. 

2.3  UNSTABLE  RESONATOR  DATA 

The  Helios  CL  II  and  CL  I  lasers  were  fully  characterized  in 
previous  studies2, 3.4,10,11,13.  por  completeness,  the  CL  II  and  CL  I 
unstable  resonator  data  are  summarized  below. 


OUTCOUPLED  POWER, WATTS 


MN0R03UR  HELIOS  CL  H  LASER 

RUN  34  MODIFIED  50%  GEOMETRIC  OUTCOUPLED 
CONFOCAL  UNSTABLE  RESONATOR 


0  4  6  12  16  20  24  28 


ITERATION  NUMBER 


Figure  3  MN0RQ3UR  convergence  parameter  PlQS$  total/pchem  and  the 

multiline  outcoupled  power  calculated  in  the  optics®,  PQUt, 

versus  iteration  number  for  rotational  nonequilibrium  for  a 
50#  geometric  outcoupled,  confocal,  unstable  strip  resonator 
with  a  large  mirror  diameter, of , 1.0  cm  for  the  Helios  CL  II 
laser  for  Run  34  at  5.4  Torr10’  .  The  FPSR  baseline 
parameters  were  used. 


OUTCOUPLED  POWER, WATTS 
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MN0R03UR  HELIOS  CL  H  LASER 

RUN  34  MODIFIED  50% GEOMETRIC  OUTCOUPLED 
CONFOCAL  UNSTALLE  RESONATOR 


0  4  8  12  16  20  24  28 


ITERATION  NUMBER 


Figure  4  MN0R03UR  outcoupled  power8,  Pout,  for  individual  lines  versus 

Iteration  number  for  a  50%  geometric  outcoupled,  confocal, 
unstable  strip  resonator  with  a  large  mirror  diameternof . 1.0 
cm  for  the  Hellos  CL  II  laser  for  Run  34  at  5.4  TorriU*  . 

The  FPSR  baseline  parameters  were  used. 


MN0R03UR  HELIOS  CL  H  LASER 

Run  34  modified  50%  geometric  outcoupled 

confocol  unstable  resonator 

NF«  9.457  Iteration  26 

OA  «  0.50  cm 


scraper  mirror 


Figure  5  MN0R03UR  gains  at  Iteration  26  for  a  50%  geometric  out- 
coupled,  confocal,  unstable  strip  resonator  with  a  large 
mirror  diameter  of. 1.0. cm  for  the  Helios  CL  II  laser  for 

Dim  A  ^  TVia  Kacol  ino  naramotorc  wc 


Run  34  at  5.4  Torr 
used. 


The  FPSR  baseline  parameters  were 
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Line 

Period 

Frequency 

Amplitude  Modulation 

Power 

ns 

MHz 

'  Pv  (J) 

Watts 

Total  Power 

40 

25 

36* 

11.00 

Pj(4) 

33 

30 

82 

1.70 

Pj(S) 

40 

25 

100 

1.20 

Pj(6) 

40 

25 

76 

1.90 

Pj(7) 

33 

30 

18 

1.00 

PjIS) 

0.05 

P2(4) 

40 

25 

27 

0.60 

P2(S) 

33 

30 

83 

2.70 

P2(6) 

40-47 

21-25 

92 

1.30 

P,(7) 

33-40 

25-30 

40 

0.50 

P 2(8)  0.05 


♦Amplitude  Modulation  for  Total  Power:  %  P^ 


Table  2  MN0R03UR  CL  II  time-dependent  oscillation  periods, 

frequencies,  and  amplitude  modulations  (%  PV(J))  for  a  50% 
geometric  outcoupled,  confocal  unstable  resonator  with  a  5  mm 
slit  for  the  Run  34  flow  rates  at  5.4  TorriU,ii.  The  FPSR 
baseline  parameters  were  used. 
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Figure  6  shows  the  layout  of  the  50%  geometric  outcoupled, 
confocal,  unstable  resonator  used  on  both  the  Helios  CL  II  and  CL  I 
lasers.  The  resonator  consists  of  a  convex  and  a  concave  mirror 
separated  by  100  cm,  with  a  flat  scraper  mirror  placed  at  an  angle  of 
45°  to  the  optical  axis  near  the  convex  mirror.  The  radius  of  curvature 
of  the  mirrors  is  -200  cm  for  the  convex  mirror  and  400  cm  for  the 


concave  mirror  which  result  in  a  resonator  magnification  of  2.  The 


scraper  mirror  consists  of  two  flat  mirrors  which  slide  apart  to  form  a 


slit  whose  width  could  be  varied  from  0.0  to  1.0  cm.  Since  the 


resonator  was  a  confocal  unstable  resonator  with  a  magnification  of  2, 


the  effective  diameter  of  the  convex  mirror  is  equal  to  the  scraper 


mirror  slit  width,  and  the  effective  diameter  of  the  concave  mirror  is 


equal  to  twice  the  scraper  mirror  slit  width.  The  resonator  Fresnel 


number  is  given  by 


*  4xL 


(2.3-1) 


where  D  is  the  effective  diameter  of  the  large  mirror,  \  is  the  optical 


wavelength  and  L  is  the  mirror  spacing.  Since  the  effective  diameter  of 


the  large  mirror  is  determined  by  the  scraper  mirror  slit  width,  the 


Fresnel  number  is  controlled  by  the  scraper  mirror  slit  width. 


2.3.1  RUN  34,  LOW  PRESSURE  DATA 


Total  power,  oscillation  frequency  and  amplitude  measurements  for 


the  low  SFg  and  low  H2  flow  rates.  Run  34,  were  taken  with  a  laser 


cavity  pressure  of  5.3  Torr  in  the  CL  II  and  5.4  Torr  in  the 
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CL  i 13^  Because  the  resonator  was  symmetric  and  was  always 
aligned  with  the  upsteam  edge  at  the  H2  injectors,  the  distance  between 
the  optical  axis  and  the  H2  injectors  was  always  equal  to  the  scraper 
mirror  slit  width.  The  oscillation  and  total  power  data  are  summarized 
as  a  function  of  Fresnel  number  for  the  CL  II  laser  in  Table  3  and  for 
the  CL  I  laser  in  Table  4.  From  Tables  3  and  4  it  is  seen  that,  in 
addition  to  a  40  nsec  oscillation,  there  was  a  7  nsec  oscillation 
superimposed  on  it.  Since  the  mirror  spacing  of  the  resonator  was  100 
cm  and  the  round  trip  transit  time  for  a  photon  is  6.67  nsec,  the  7  nsec 
oscillation  corresponds  to  a  mode  beat  of  the  laser.  The  time-dependent 
oscillations  did  not  occur  for  Fresnel  numbers  below  a  demarcation 
Fresnel  number  of  1.5  in  the  CL  II  and  3.2  in  the  CL  I.  Since 
calculations  had  shown  that  the  length  of  the  saturated  gain  region  was 
independent  of  the  size  of  the  resonator-*,  Tables  3  and  4  show  that  the 
amplitude  of  the  oscillations  increased  as  the  fraction  of  the  resonator 
filled  by  the  oscillating  lines  decreased. 

The  oscillation  data  for  individual  lines  as  a  function  of  Fresnel 
number  are  given  for  the  CL  II  and  CL  I  in  Tables  5  and  6, 
respectively.  These  data  show  that  the  40  nsec  oscillation  on  lines 
whose  saturated  gain  did  not  fill  the  resonator  did  not  occur  for 
Fresnel  numbers  below  the  demarcation  Fresnel  number  of  1.5  for  the  CL 
II  and  3.2  for  the  CL  I. 


2.3.2  RUN  36,  LOW  PRESSURE  DATA 


Oscillation  frequency  and  amplitude  measurements  for  Run  36,  the 
middle  SFg  and  high  H2  flow  rate,  were  taken  at  laser  cavity  pressures 
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CL  II,  Run  34,  5.30  Torr,  50%  Geometric  Outcoupled, 
Symmetric,  Confocal,  Unstable  Resonator 


Scraper  Mi rror  Nc 
Slit,  am 


p 

rcav 

(Torr) 


PT  OSC. 

Period 
(Watts)  (nsec) 


Amplitude 

Modulation 


1.0 

0.357 

5.31 

7.25 

7 

2.0 

1.428 

5.30 

11.7 

7 

— 

3.0 

3.214 

5.29 

9.5 

40/7 

3.0 

4.0 

5.714 

5.31 

7.0 

40/7 

22.0 

5.0 

8.929 

5.34 

2.25 

40/7 

25.0 

6.0 

12.857 

5.31 

0.98 

40 

50.0 

Table  3  Measured*®’**  frequency  and  amplitude  of  the  time-dependent 
oscillation  of  total  power  as  a  function  of  Fresnel  number 
for  the  Run  34  flow  rates  at  5.3  Torr. 
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CL  I,  Run  34,  5.40  Torr,  50X  Geometric  Outcoupled 
Symmetric,  Confocal,  Unstable  Resonator 


Scraper  Mirror  NF  PT  OSC.  Amplitude 

Slit,  m  Period  Modulation 

(Matts)  (nsec) _ (X  Py) _ 


1.0 

0.357 

0.39 

2.0 

1.428 

1.97 

7 

5.3 

3.0 

3.214 

1.03 

7 

4.7 

4.0 

5.714 

0.10 

47/7 

50.1 

5.0 

8.929 

0.01 

Table  4  Measured^  frequency  and  amplitude  of  the  time-dependent 

oscillation  of  total  power  as  a  function  of  Fresnel  number 
for  the  Run  34  flow  rates  at  5.4  Torr. 


Cl  II,  fcun  34  S0\  Geometric  Out  coupled,  S>Tu»etnc,  Confoc.il,  Unstable  Resonator 

Cavity  Pressure  3.30  Torr 


>o 

n 

a 


r*  rs 


TJ  N  C 

O  U  o 

'  '  T)  - 
3  *- 
«-■  « 


Ol 

o» 


O  Ci  ul 
V  l«N 


e*  © 
^  <« 

O  Cl  • 

«e  a  *o 


»■>-*© 

O  t  . 

^  a  »*» 


w  a 

o  o 


£  3 

«• 

o  o 


*  o 

if 

a.  x 

o  ta 
to  c 
O  4» 


—  3 
0.-0  - 

ii 


u 

l/i 

O 


c  C  X 
■  oX 

T»  W  * 
O 

~  <r 

»-  o  c 
V  u 

a  x  u. 

u  u  J 
lO  C  tO 

O  v  O 


o 

•o 


O  i 

^  o 


O  <9  I 

<«r  &>  o 
a  r* 


o 

*o 


o 

to 

o 


c  “O 

o  ^  o  »*« 

—  •«*  o 


3  w  w 
*-  *  > 
-  «  a. 
— >  o 
a.  T3  ^ 

5£ 


u  C  u 
t#  1/) 

O  u  o 


24 


CL 

5 


—  a 

3 

T3  ** 

O 


U  •*-  C  T3  •—  O 


c  > 

o 

01  <d 

01 

"U 

4-> 

a;  •»- 

•r— 

C  3 

4-> 

01 

3  XJ 

4-> 

—  X3 

Ol 

an  <u 

CT  C 

U 

*3  «r- 

X3 

(O 

*-> 

ai 

c 

4-*  > 

s- 

ITS 

v. 

3 

-Q  *r- 

4-5 

01 

V- 

4-  C 

4— 

o  -o 

00 

> 

aj 

O 

c 

to 

to 

CL  • 

-H 

<d 

4-»  *r- 

4- 

o  s- 

— «  00 

o 

01 

o 

•*  c 

00 

C  01 

01 

E 

O  4-> 

D  O 

<d 

JC 

JZ 

4->  O 

-4  •»— 

01  4-5 

4-» 

4-» 

ai 

“O  4-» 

u 

3- 

j«:  4-* 

01  CT3 

u 

01  01 

01 

01 

in  <u 

V-  i— 

o 

5  VI 

4-> 

s- 

a»  xj 

3  e— 

H- 

3 

ra 

3 

3 

l/>  *e- 

»d  <d 

V- 

oo 

qj 

<d  U 

ro 

4-»  U 

01 

*d 

o  1/1 

01  00 

• 

to  01 

a* 

ai 

O  -O 

x  o 

ir> 

X>  -Q 

o 

E  U  Ol 

cr> 

01 


-a 

*3 


<o 

o 

ai 

Ss 

4-> 

4-5 

JC 

J1.’ 

C  L0 

01 

00  4-5 

QJ  (U 

c 

J* 

01 

-o  +-> 

•r* 

Ol 

fd 

C  01 

C  i« 

r— 

_d 

01 

O  -r-  4-5 

ou  s- 

■M 

4->  -J  <0 

Xy 

Q. 

r— 

-o 

U 

O 

ai  » 

tO 

■r— 

o 

TD  01 

•-* 

-o  o 

3 

3E 

o 

01  •  CL 

u 

4> 

i  r- 

XJ 

4-» 

00  V)  o 

*-* 

QJ  4- 

•r— 

4-> 

3  C 

a> 

E 

> 

*f~ 

01 

O  XI 

Vj 

•r-  KT 

•r— 

i — 

3- 

00  •»-  r— 

sH 

p 

4J  ro 

-o 

1/1 

01 

to  4->  3 

r J 

v> 

c 

S 

£30 

«•* 

QJ  c 

>— ♦ 

E 

J3  £ 

JC  3 

E 

CO 

3-  t- 

rr 

a> 

4->  Q£ 

01 

O  W  3- 

*-* 

• 

v£> 

oo 

4-»  4-5  O 

yj 

4-  01 

u 

<d 

(J  CO  4-» 

Ml 

o 

O  JC. 

01 

-o 

O 

01  *r-  O 

<« 

+J 

-Q 

c 

4->  d  01 

QJ 

E 

hj 

*01 

01  4-5 

Cl 

n. 

TO  S- 

3 

00 

"TD  t—  01 

■j,' 

O 

=J  O 

C 

E 

01 

<d  x? 

>3 

4J  4- 

E 

JC 

01  s- 

O 

•r- 

4-5 

W)  4-5  CO 

f—  1/1 

01 

in 

o  u< 

J4 

CL  QJ 

c 

c 

CL  QJ  C 

f8 

c 

E  C 

CO 

aj 

CL*—* 

*1 

iO  ’i- 

01 

JZ 

c  CO 

r*“ 

3- 

4—1 

VI 

4-5 

o 

o 

-a 

U_ 

01 

3-  oo 

w 

C  i— 

I- 

c 

•  01  tO 

<o  <a 

4- 

o 

•r— 

3-  £  4- 

a 

3 

O 

4- 

O  O 

a 

>lTJ 

4-5  a.  oi 

S:  ' 

.T*  j 


K  ' 

vi  : 


§ 


3  i 

S  1 


s  I 


•5*  < 
■y  ' 

-V  » 


o 

u 

i  I 


" 


•*  - 
a  — 
m  oi 


* 

V 


S'  '-‘.'tvW 


w-KT-jr.! 


Run  34,  5.40  Torr,  50%  Geometric  Outcoupled, 
Symmetric,  Confocal,  Unstable  Resonator 
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of  6.30  Torr  In  the  CL  II  and  6.50  Torr  in  the  CL  il0.ll.13.  Since  a 
symmetric  resonator  that  was  aligned  with  the  upstream  edge  at  the  H2 
injectors  was  used  in  both  cases,  the  distance  between  the  optical  axis 
and  the  H2  injectors  was  always  equal  to  the  scraper  mirror  slit 


width.  These  data  are  summarized  in  Tables  7  through  10.  Again,  both 
the  7  nsec  and  the  40  nsec  oscillations  appear  in  both  lasers  on  lines 
whose  saturated  gain  does  not  fill  the  resonator.  These  oscillations 


did  not  occur  for  Fresnel  numbers  below  the  demarcation  Fresnel  numbers 
of  1.5  for  the  CL  II  and  3.2  for  the  CL  I. 


2.3.3  COMPARISON  OF  CL  II  AND  CL  I  UNSTABLE  RESONATOR  DATA 

It  should  be  noted  that  a  comparison  of  the  CL  II  and  CL  I  unstable 
resonator  data  showed  that  fewer  lines  were  observed  with  the  CL  I  and 
the  total  powers  of  the  two  devices  are  very  different.  In  Ref.  13, 
this  is  shown  to  be  a  consequence  of  the  fact  that  the  saturated  gain  in 
the  CL  I  is  twice  as  large  as  the  saturated  gain  in  the  CL  II.  Overall, 
as  shown  in  Ref.  13,  the  CL  I  unstable  resonator  data  are  consistent 
with  the  CL  II  unstable  resonator  data. 


2.4  COMPARISON  OF  THE  A  PRIORI  UNSTABLE  RESONATOR  CALCULATIONS  TO 

THE  CL  II  DATA 

Tables  2,  3  and  5  show  that  when  compared  to  the  data,  the  5  mm 
slit  MN0R03UR  calculation  in  Section  2.2  accurately  predicted  both  the 
period  and  amplitude  modulation  of  the  40  nsec  time-dependent  oscilla¬ 
tions.  However,  the  total  power  was  overpredicted  beyond  what  was 
expected  from  the  lack  of  Brewster  window  losses  in  the  calculation. 
This  Indicated  that  the  FPSR  baseline  parameters  used  in  the  calculation 
were  not  correct. 
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CL  II,  Run  36,  6.30  Torr,  50%  Geometric  Outcoupled, 
Symmetric,  Confocal,  Unstable  Resonator 


Scraper  Mirror  NF 
Slit,  mm 


pcav 

(Torxl 


PT  OSC. 

Period 
(Watts)  (nsec) 


Amplitude 

Modulation 


1? LPt1„ 


1.0 

0.357 

6.5 

10.5 

-- 

2.0 

1.428 

6.5 

17.2 

7 

— 

3.0 

3.214 

6.5 

12.2 

40/7 

3.7 

4.0 

5.714 

6.5 

6.55 

40/7 

33.0 

5.0 

8.929 

6.5 

1.55 

40/7 

37.5 

Table  7  Measured*®’^  frequency  and  amplitude  of  the  time-dependent 
oscillation  of  total  power  as  a  function  of  Fresnel  number 
for  the  Run  36  flow  rates  at  6.3  Torr. 
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CL  I,  Run  36,  6.50  Torr,  50%  Geometric  Outcoupled, 
Symmetric,  Confocal,  Unstable  Resonator 


Scraper  Mirror  NF 
Slit,  mm 


(Watts] 


OSC.  Amplitude 
Period  Modulation 
(nsec)  (X  PT) 


1.0 

0.357 

1.08 

— 

-- 

2.0 

1.428 

2.59 

7 

3.7 

3.0 

3.214 

1.22 

7 

2.9 

4.0 

5.714 

0.11 

47/7 

15.5 

5.0 

8.929 

0.01 

40 

80.8 

1  1 

Table  8  Measured1-3  frequency  and  amplitude  of  the  time-dependent 

oscillation  of  total  power  as  a  function  of  Fresnel  number 
for  the  Run  36  flow  rates  at  6.5  Torr. 
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Run  36,  6.50  Torr,  50%  Geometric  Outcoupled, 
Symmetric,  Confocal,  Unstable  Resonator 


In  addition  to  the  40  nsec  oscillation,  the  data  showed  a  7  nsec 
oscillation  superimposed  on  it.  Since  the  mirror  spacing  of  the 
resonator  was  100  cm  and  the  round  trip  transit  time  for  a  photon  is 
6.67  nsec,  the  7  nsec  oscillation  corresponds  to  a  longitudinal  mode 
beat  of  the  laser.  The  MN0R03UR  calculation  did  not  predict  a  7  nsec 
oscillation  because  the  model  does  not  include  the  longitudinal  modes  of 
the  lasing  lines. 

2.5  UNSTABLE  RESONATOR  BASELINING  CALCULATIONS 

The  first  step  required  to  determine  if  it  was  necessary  to 
baseline  the  unstable  resonator  model  to  unstable  resonator  data  was  to 
use  the  Fabry-Perot  and  stable  resonator  parameters  in  the  unstable 
resonator  model.  If  calculations  with  these  parameters  did  not 
correctly  model  the  measured  unstable  resonator  performance,  the 
unstable  resonator  model  would  need  to  be  baselined  to  unstable 
resonator  data.  Unstable  resonator  calculations  with  the  Fabry-Perot 
and  stable  resonator  parameters  were  performed  for  two  50%  geometric 
outcoupled,  confocal  unstable  resonators,  one  with  a  4  mm  slit  and  one 
with  a  5  mm  slit.  The  resonator  with  the  4  mm  slit  was  the  resonator 
for  which  the  best  CL  II  time-dependent  oscillation  data  were 
available*®’**.  The  resonator  with  the  5  mm  slit  was  the  one  for  which 
the  a  priori  calculation  was  performed  that  predicted  the  existance  of 
time-dependent  oscillations  on  lines  whose  saturated  gain  does  not  fill 
the  resonator. 

The  results  of  these  calculations  are  summarized  in  Table  11. 
Results  from  calculations  with  the  FPSR  parameters  are  also  included  in 


Baseline 

Parameters 


Slit 

Width 


Baseline 

Parameters 

Slit 

Width 

(ran) 

Total 

Power 

(Watts) 

Oscillation 

Period 

(nsec) 

Oscillation 
Amplitude 
(%  total  power) 

S 

EMM  Data10’11 

4 

7.0 

40 

22.0 

s? 

y 

_  •  f 

; 

Fabry-Perot 

4 

19.9 

20 

2.5 

IL 

Stable  Resonator 

4 

30.4 

20 

1.6 

i\ 

‘ 

FPSR 

4 

11.6 

20 

1.7 

£  : 

EMM  Data10’11 

5 

2.25 

40 

25.0 

§ 

• 

Stable  Resonator 

5 

29.0 

20 

1.7 

FPSR 

5 

11.0 

40 

36.0 

£  i 

*  m 

Table  11  Summary  of  the 
MN0R03UR  calcul 
pressure  Run  34 

4  mm  and  5  mm 
ations  for  the 
flow  rates. 

slit  unstable  resonator  data  and 
Helios  CL  II  for  the  low 

All  calculations  used  tw  =  1.0. 
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Table  11  since  these  parameters  were  used  in  previous  unstable  resonator 
modeling.  All  of  the  calculations  summarized  in  Table  11  had  no 
Brewster  window  losses  included  and  thus  were  expected  to  overpredict 
the  total  power.  Table  11  shows  that  the  FPSR  parameters  allow  the 
unstable  resonator  model  to  more  accurately  predict  the  data  than  either 
the  Fabry-Perot  or  stable  resonator  parameters.  However,  with  the  FPSR 
parameters,  the  model  overpredicted  the  total  power  in  both  resonators 
beyond  what  was  expected  from  the  lack  of  Brewster  window  losses  in  the 
calculation.  Also,  the  period  and  amplitude  modulation  of  the  time- 
dependent  oscillations  were  not  correct  in  the  4  mm  slit  case.  Based 
upon  the  fact  that  neither  the  Fabry-Perot,  nor  the  stable  resonator, 
nor  the  FPSR  parameters  allowed  the  unstable  resonator  model  to 
accurately  predict  unstable  resonator  performance  as  the  size  of  the 
resonator  was  varied,  the  unstable  resonator  model  was  baselined  to  the 
unstable  resonator  data. 

To  baseline  the  unstable  resonator  model,  four  baseline  parameters 
had  to  be  determined  so  that  the  unstable  resonator  model  tracked  the 
unstable  resonator  data  as  resonator  size,  flow  rates  and  size  of  the 
laser  were  varied.  The  data  to  be  tracked  were  the  total  power  and  the 
period  and  amplitude  modulation  of  the  time-dependent  oscillations  that 
occur  on  lines  whose  saturated  gain  does  not  fill  the  resonator.  The 
four  baseline  parameters  are  the  static  temperature  of  the  flow  at  the 
H2  injectors,  the  percent  dissociation  of  the  SFg  in  the  arc,  the 
primary  mixing  length,  and  the  secondary  mixing  length.  Of  these  four 
parameters,  two  had  been  fixed  with  confidence  in  previous  studies* 

A  primary  mixing  length  of  0.001  cm  and  a  temperature  of  the  flow  at  the 


Hg  Injectors  of  450  K  were  found  to  be  appropriate.  Thus,  only  two  of 
the  four  parameters,  the  percent  dissociation  of  the  SFg  in  the  arc  and 
the  secondary  mixing  length,  were  varied  in  baselining  the  unstable 
resonator  model. 

In  an  earlier  investigation*,  three  of  the  four  baseline  parameters 
were  varied  and  the  effect  on  the  total  power,  power  split  between 
vibrational  bands  and  beam  diameter  was  calculated  using  the  rotational 
equilibrium,  Fabry-Perot  BLAZE  II  model.  The  one  baseline  parameter 
that  was  not  varied  in  this  parametric  study  was  the  primary  mixing 
length  which  was  set  at  0.001  cm.  The  results  of  these  calculations, 
which  served  as  a  guide  in  choosing  parameter  combinations  for  the 
unstable  resonator,  are  presented  in  Figs.  7  through  12. 

The  low  pressure  Run  34  flow  rates  were  chosen  as  the  baseline  flow 
rates  since  the  best  unstable  resonator  data  were  taken  at  these  flow 
rates*®’**.  The  2  mm  slit  resonator  was  chosen  as  the  first  baseline 
case  because  experimentally  it  does  not  display  the  time-dependent 
oscillations.  Once  the  data  were  matched  for  the  2  mm  slit,  the  data 
for  larger  slit  sizes  would  be  matched  if  the  parameter  set  was  valid. 
From  Figs.  7  through  12,  a  set  of  parameters  was  chosen  to  match  the  2 
mm  slit  total  power  data.  This  parameter  set  is  denoted  Set  A  in 
Table  12.  Two  calculations  were  made  with  Set  A,  one  with  Brewster 
window  losses  included  and  one  with  no  Brewster  window  losses.  This 
provided  a  check  on  the  validity  of  using  the  stable  resonator  effective 
Brewster  window  transmissivity  of  tw  =  0.97  in  the  unstable  resonator 
model.  Table  13  summarizes  the  results  of  the  2  mm  slit  calculations 
with  Set  A.  A  comparison  of  the  predicted  total  powers  indicated  that 
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HELIOS  CL  IT  LASER 
BLAZE  IE  r  •  0.57 
Run  36  Low  Pressure 
Lmp«  0.001cm 

■  %  SF6  diss*  4.00iT,«  450*K 
▲  %SF6di*s»  2.75,  T(«  450* K 
•  %  SF6  diss  «  2  00,  T,  «  450*  K 
T  %  SF6  diss  «  2.75-,  T,*550*K 


Figure  8 


Variation  of  lasing  cutoff1  as  a  function  of  secondary  mixing 
length  from  the  8LAZE  II  model  of  the  CL  II  laser.  L^p  is 
the  mixing  length  of  the  primary  stream. 
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HELIOS  CL H  LASER 
BLAZE  H  r  ■  0.57 
Run  36  Low  Pressure 


%SF6  Dissociotion 


Figure  9  Variation  of  total  power  and  fraction  of  power  in  the  1*0 
vibrational  band1  as  a  function  of  precent  SFg  dissociation 
from  the  BLAZE  II  model  of  the  CL  II  laser,  LMp  is  the  mixing 
length  of  the  primary  stream,  L^j  Is  the  mixing  length  of  the 
secondary  stream. 


HELIOS  CL  H  LASER 
BLAZE H  r  *  0.57 
Run  36  Low  Pressure 
Lmp»  0.00 1  cm 

▼  Lms«  7  5cm  jT,  «  450*K 
#  LMs  *  6.5cm;  Tj  *  450#K 
ALMS-SScmjT,  *  450*K 
■  Lms * 9.5 cn^T,  •  350*K 


%SF6  Dissociotion 


Figure  10  Variation  of  lasing  cutoff  as  a  function  of  percent  SFg 

dissociation  from  the  BLAZE  II  model  of  the  CL  II  laser.  Lj, 
Is  the  mixing  length  of  the  primary  stream.  LM$  is  the 
mixing  length  of  the  secondary  stream. 


Total  Power  (watts) 


Tj  ot  H2  Injectors  (°K) 

Figure  11  Variation  of  total  power  and  fraction  of  power  in  the  l-»0 

vibrational  band1  as  a  function  of  the  initial  temperature  at 
the  H2  Injectors  from  the  BLAZE  II  model  of  the  CL  II 
laser.  L^n  Is  the  mixing  length  of  the  primary  stream.  Luc 
Is  the  mixing  length  of  the  secondary  stream. 
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Figure  12  Variation  of  the  lasing  cutoff1  as  a  function  of  the  initial 
temperature  at  the  H2  injectors  from  the  BLAZE  II  model  of 
the  CL  II  laser.  L^p  is  the  mixing  length  of  the  primary 
stream,  L^<j  is  the  mixing  length  of  the  secondary  stream. 
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Parameter  Set  Static  Temp,  at  Primary  Secondary  X  SFg 

Ho  Injectors  Mixing  Mixing  Dissociation 

(K)  Length  Length 

(cm)  (cm) 


Baseline  Slit  Total  Oscillation  Oscillation 

Parameters  Width  Power  Period  Amplitude 


Set  A  (tw  =  1.0) 

2 

26.2 

20 

7 

•v 

Set  A 

2 

19.6 

20 

7 

3 

Set  B 

2 

8.8 

20 

8 

£ 

Set  C 

2 

12.0 

20 

8 

Set  C 

4 

20.4 

40 

9 

% 

Set  0 

2 

4.7 

-- 

— 

Set  D 

4 

7.3 

40 

9 

'? 

Set  E 

2 

5.3 

40 

4 

■ 

Set  E 

4 

7.5 

40 

4 

y 

Table  13  Summary  of  2  mm  and  4  mm  slit  MN0R03UR  baselining  calculations  $ 
and  data  for  the  Helios  CL  II  for  the  low  pressure  Run  34  flow 
rates.  All  calculations  were  with  tw  =  0.97  unless  otherwise 
Indicated;  all  calculations  were  for  a  50%  geometric  outcoupled,  ft 
confocal,  unstable  strip  resonator. 


tw  *  0.97  Introduced  a  25%  loss  which  compared  well  with  the  27%  loss 
Introduced  in  the  stable  resonator  model.  Table  1.  The  comparable 
reduction  in  total  power  in  both  stable  and  unstable  resonator  models 
reinforced  the  conclusion  that  the  value  of  the  effective  transmissivity 
of  the  Brewster  windows  determined  for  the  stable  resonator  model  is 
appropriate  for  the  unstable  resonator  model. 

Table  13  shows  that  the  calculation  with  Set  A  that  included 
Brewster  window  losses  overpredicted  the  total  power.  Since  Figs.  9  and 
10  Indicated  that  total  power  is  a  strong  function  of  the  percent 
dissociation  of  SFg  in  the  arc,  to  reduce  the  predicted  total  power,  the 
SFg  dissociation  was  reduced  to  2.36%.  This  new  parameter  set  is 
denoted  Set  B  in  Table  12.  The  results  of  a  2  mm  slit  calculation 
performed  with  Set  B  are  shown  in  Table  13.  The  total  power  calculated 
with  Set  B  was  too  low. 

A  value  of  the  percent  SFg  dissociation  that  would  give  the 
observed  2  mm  slit  total  power  was  determined  through  linear  inter¬ 
polation  from  the  calculation  with  Set  A  and  the  calculation  with 
Set  B.  This  percent  SFg  dissociation,  3.2%,  along  with  the  other 
parameters  from  Set  B,  comprise  parameter  Set  C,  Table  12.  A  2  mm 
unstable  resonator  calculation  with  parameter  Set  C,  Table  13,  gave  good 
agreement  with  the  measured  total  power.  There  was,  however,  a  20  nsec 
small  amplitude  oscillation  predicted  that  was  not  seen  in  the  data. 

Since  2  mm  slit  calculations  with  parameter  Set  C  matched  the  total 
power  data,  a  calculation  with  the  same  parameter  set  was  made  for  a  4 
mm  slit  resonator.  Table  13.  Unfortunately,  with  parameter  Set  C,  the 
4  mm  slit  total  power  was  overpredicted  by  a  factor  of  three.  The 


oscillation  period,  however,  was  correct  at  40  nsec.  The  inability  of 
Set  C  to  allow  the  model  to  predict  the  data  as  the  resonator  size 
changed  indicated  that  Set  C  was  not  correct. 

The  gains  and  the  Intensity  distributions  from  the  preceding  two 
calculations  with  Set  C,  Figs.  13  through  16,  were  compared.  These 
figures  suggested  that  the  gain  zones  should  be  shortened.  With  the 
gain  zones  shortened,  the  intensity  peak  would  move  out  from  behind  the 
small  mirror,  making  It  less  symmetric  in  the  2  mm  slit  case. 
Laboratory  burn  blocks  show  an  unsymmetric  intensity  distribution  in  the 
2  mm  slit  case.  Since  moving  the  intensity  peak  into  the  outcoupling 
portion  of  the  cavity  would  Increase  the  total  power,  the  percent  SFg 
dissociation  should  be  reduced  to  compensate. 

To  choose  a  parameter  set  that  would  shorten  the  gain  zones,  a  set 
of  BLAZE  II  zero  power  gain  calculations  were  performed  to  determine  the 
variation  of  the  gain  distribution  as  a  function  of  the  secondary  mixing 
length  and  the  percent  SFg  dissociation.  The  resulting  zero  power  gain 
curves  are  presented  in  Figs.  17  through  30.  The  BLAZE  II  zero  power 
gain  curves  indicated  that  the  magnitude  of  the  peak  zero  power  gain  was 
a  strong  function  of  the  percent  SFg  dissociation  and  a  very  weak 
function  of  the  secondary  mixing  length  and  that  the  axial  location  of 
the  peak  zero  power  gain  is  a  strong  function  of  the  secondary  mixing 
length.  Figures  29  and  30  indicated  that  a  secondary  mixing  length  of 
2.5  cm  would  move  the  gains  upstream  and  a  1.2%  SFg  dissociation  would 
give  peak  zero  power  gains  of  about  0.08  cm"1.  This  peak  zero  power 
gain  and  position  were  reasonable  and  thus  another  set  of  baseline 
parameters  was  chosen.  These  parameters  are  denoted  Set  D  in  Table  12. 
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Figure  23  Variation  of  the  zero  power  gain  as  a  function  of  distance 
downstream  from  the  injectors  from  the  BLAZE  II  model  of 
the  CL  II  laser.  L  is  the  mixing  length  of  the  primary 
stream,  L  is  the  mixing  length  of  the  secondary  strr  i,  T 
is  the  static  temperature  of  the  flow  at  the  H?  injectors. 


0.20  | - - 

m 

0.18-  Helios  CLII  Laser,  Blaze  II 


c  41  x  sz  a> 

•f-  I_  -r-  4->  -C 

<0  O  E  CT14-J 

CT14J  C 

0  01014- 

i-  aj  -c  >—  o 

Ol  1->4-» 

5  C  CD  OJ 

O  41  C  t- 

cx  •>—  *i—  rj 

<\j  X  4-> 

o  x  Oj-  <o 
u  EE  1- 
Ol  Ol  — I  O) 
(M  -C  Ol  Q. 

+->  X  E 

ai  .  4->  oi 

•C  E  s-  4-1 
4-1  O  Ol  I/I 

S-  41  -i-  O 
4-  4-  id  ■<- 

o  i —  l/4-i 

E  E<d 

C  rQ  •—  _ I  4J 

0  0)“  41 

•■—l-  - 

4+1  J  E  01 
«  4iu 
•—  C  Ol  4-1 

t  J  01  l. 

Id  O  X  4->  41 

>  "O  4-»  41  •<- 


CD 

OJ 

o 

co 

o 

CD 

o 

q 

0J 

o 

o 

o 

O 

O 

o 

o’ 

o 

o' 

o 

o 

d 

((_.  LUO)  UIBQ 


v.v.v.v.v;  /  ’ -v -v  ^  *:•. 


rc» 


Figure  26  Variation  of  the  zero  power  gain  as  a  function  of  distance 
downstream  from  the  H2  injectors  from  the  BLAZE  II  model  of 
the  CL  II  laser.  L  is  the  mixing  length  of  the  primary 
stream,  L  is  the  mixing  length  of  the  secondary  stream,  T 
is  the  static  temperature  of  the  flow  at  the  Ho  injectors. 


,1.  tf 


5  <1)  r  E  o 

3  ^  0  O  O 

w  £  lo  cm  o 

rT  m  P  cO  co  u 

5  ^  Q.  II  II  c 


CO  v  ID  CD  N 


”3  “3  “3  “3  “3  “ 


****** 


=  O  N 


CO  O 

E  ° 

*—  //> 


«  Q 


s  §^3 

§  <*  "  CO  ■ 
CC  Q_  >  H- 


t*- 

1— 

CD 

o 

• 

U 

>> 

•»  CO 

c 

r—  S- 

E  U 

<d 

<u  <a 

«3  o 

4-> 

•o  E 

at  -t-4 

c/> 

o  -r- 

u  u 

•f— 

E  J- 

4->  tu 

T3 

o.  co  *n 

►— « 

c 

M- 

►— *  dj 

o 

_c 

u 

LlJ  4-> 

Id  CM 

c 

r-M 

-a  =n 

O  <C  «- 

c 

•r— 

— 1  o 

O  CD 

4— >  CO 

U  -C 

o 

-C 

CD  4-> 

c 

CD  -♦-> 

CO 

3 

sz  cn 

4-> 

M— 

4— >  C 

CD  Id 

HK9 

-C 

1 _ 1  ^ 

«u 

o 

4-»  JC 

o 

to  U 

0*4-  r— 

03  4- 

c  o  u- 

•r— 

C  to 

X  -C  <D 

V- 

•*-  4->  -C 

>o  o 

E  CT4-> 

C 

u 

<D  <D  M_ 

V-  CL) 

-C  i—  o 

CL)  T-> 

4-> 

3;  C 

cn  <d 

il  O  (U  U1 

i-  to  f-  (J 

U-  U—  ra  ••— 
O  » —  i/u-> 

E  E<o 
c  <d  •— i  — i  *j 
o  <L)  •— I  1/1 
•f—  V-  » 

4->  4->  _J  E  <U 

<C  to  c_>  n3 

•i-  C  oi  -u 

<e  o  .c  4->  to 

>  T3  -M  to  -r- 


o 

CM 

00 

CO 

CM 

o 

co 

o 

CO 

o 

o 

CM 

o 

o 

o 

o 

d 

o 

o 

o 

d 

o 

o 

d 

o 

o 

(  LUO)  UjEQ 


co  ^  lo  co  oo 

ii  ii  ii  ii  ii  ii 

“D  "3  “3  "0  “D  “ D 


0  I  +  M  ■!» 


U  >,  »  </> 

C  .—  t.  E  l_ 

«  ill  n)  «  o 

+JTJ  E  ai*> 
H  O  'i-  L  U 

•i-  E  i-  +->  CL) 

"O  Q-  l/l  •'—> 


C/3 

—  c 


4-  •— I  CD  >,•!- 

o  .e  S- 

LlJ  -w  <d  cm 
CN  T3I 

O  <£  4-  C 
•I-  — I  O  O  CD 
4->  CD  ur 

U  JZ  (D  4-> 
C  QJ  4->  on 
D  r  oi  -m 
4-  4->  C  CLI  id 
CD  JC 


CD  CD 

N  O 

CD  CD  ^ 

m  □  0 
.  w  =S 
1-  c/3  o 

CD  CD  Q_ 
CD  i- 
CD  CL  O 

-1  03 

=  O  N 


§  8 

LO  C\J  o 

CM  CM  n 

II  II  c 

C/3  o  E 

E  o  _j 


0-15 

o  J  E 


-Q  ^ 
o  o 
t  coin 
I  LL 

T  co  II 


I  CC  Q_  >  h- 


id  Ei- 

4->  3 

O 

O 

</>  U  034-  » — 

c 

O  4- 

•r— 

C  W  X 

JZ  03 

•f-  l-  »r— 

4->  JZ 

IO  O  E 

C7>  4— > 

0>4-> 

C 

<J  03 

03  4- 

u  <v  sz 

*—  O 

0) 

X  C 

03  a; 

o  f  </> 

c  u 

Q.  •<- 

•*-  13 

CM 

X  4-> 

O  re  Q|-  Id 

S-  EE  D 

QJ  CD  — 1 

03 

N  -C 

03  Q- 

4— > 

-C  E 

<u 

4-»  03 

-C  E  s- 

4-> 

■M  O  (U 

to 

V-  1/3 

•r-  U 

*4-  4-  fd 

•r— 

O  r— 

</U-> 

E 

E^ 

C  (Om 

— J  4— > 

o  03  •— » 

(/> 

•*-  V- 

•» 

4->  4->  _J 

E  03 

<T3  CO  (_> 

«T3  -C 

•<-  c 

03  4-3 

1-  3  03 

J- 

ffl  0£ 

4->  CO 

>  “O  4-> 

CO  •*-* 

o 

00 

CD 

tj- 

CM 

o 

CO 

CD 

CM 

o 

3 

CM 

T— 

i— 

o 

o 

o 

o 

o 

•r~ 

lx. 

d 

d 

o 

o 

o 

o 

o 

o 

o 

o 

o 

(k  LUO)  UjBQ 


o.v 


•V  »'■ 


CO  n-  LO  CD  N  CO 

II  II  II  II  II  II 

$  It  H  ^ 


to 

—  c 

<d  'cd 

N  O 

Jo,. 

DQ  □  O 

i^TcO  O 

CD  CD  Q_ 

Eg  £  o 

— 1  5  CD 

=  o  N 

-1  — 1  r- 

O  -  O 

</>  S  S 

o”2 

(D  D  ^ 

Z  H  CL 


E  o 
o  o 

in  c> 

c\i  ’r“  ii 

II  II  Q_ 

W  o  E 

W  — J 

— 1  C/5 

-S* 

T  CD  S 

cli  t 


^  CO  11 

II 

vO - 

>0^-1 — 


o 

oo 

CD 

eg 

o 

CO 

CD 

eg 

o 

vC1 

CM 

”i— 

T— 

T— 

T— 

-r- 

o 

o 

o 

o 

o 

o 

d 

d 

o 

o 

o 

o 

o 

o 

o 

o 

(t.  LUO)  U|B0 


Figure  29  Variation  of  the  zero  power  gain  as  a  function  of  distance 
downstream  from  the  H2  injectors  from  the  BLAZE  II  model  of 
the  CL  II  laser.  L  is  the  mixing  length  of  the  primary 
stream,  L-,.  is  the  mixing  length  of  the  secondary  stream,  T 
is  the  static  temperature  of  the  flow  at  the  H2  injectors. 
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A  2  mm  slit  unstable  resonator  calculation  with  parameter  Set  D  was 
performed.  As  shown  in  Table  13,  the  total  power  predicted  by  this 
calculation  was  too  low  by  a  factor  of  two.  Figure  31  shows  that  the 
intensity  peak  did  not  move  as  desired.  A  comparison  of  the  2  mm  slit 
calculations  with  Sets  C  and  0  showed  that  Set  0  had  cut  the  total  power 
by  approximately  60%.  A  60%  decrease  in  the  total  power  calculated  with 
Set  B  in  the  4  mm  slit  case  would  result  in  a  total  power  of  about 
8  Watts,  which  would  agree  with  the  data.  Thus  a  4  mm  slit  calculation 
with  Set  0  was  performed.  The  results  of  this  calculation.  Table  13, 
show  that  the  4  mm  slit  data  was  matched  well  in  both  total  power  and 
the  period  of  the  time-dependent  oscillations.  The  calculated  amplitude 
modulation,  however,  was  low  compared  to  the  data. 

A  plot  of  total  power  vs.  percent  SFg  dissociation  was  constructed 
from  the  existing  calculations.  Figure  32  indicated  that  a  near  linear 
relationship  existed  between  total  power  and  percent  SFg  dissociation, 
independent  of  the  secondary  mixing  length.  This  relationship  held  for 
both  the  2  mm  and  4  mm  slit  cases.  Figure  32  also  suggested  that  if  the 
slope  of  the  2  mm  slit  line  could  be  made  larger  than  that  of  the  4  mm 
slit  line,  a  set  of  parameters  that  would  accurately  predict  both  2  mm 
and  4  mm  slit  performance  could  be  determined.  Thus,  in  an  attempt  to 
increase  the  slope  of  the  2  mm  slit  line,  the  secondary  mixing  length 
was  drastically  cut  to  0.95  cm.  The  resulting  parameter  set,  denoted 
Set  E,  is  presented  in  Table  12.  The  results  of  a  2  mm  slit  calculation 
with  parameter  Set  E  are  summarized  in  Table  13  and  Fig.  32.  This 
calculation  indicated  that  the  slope  of  the  2  mm  slit  line  did  not 
become  greater  than  that  of  the  4  mm  slit  line. 


MN0R03UR  HELIOS 
CLH  LASER 


Figure  31  MN0R03UR  total  intensity  at  iteration  18  for  a  50%  geometric 
outcoupled  confocal  unstable  strip  resonator  with  a  large 
mirror  diameter  of  4  mm  for  the  Helios  CL  II  for  the  low 
pressure  Run  34  flow  rates. 


A  final  4  mm  slit  calculation  was  made  to  investigate  whether  a 
drastic  cut  in  secondary  mixing  length  would  make  the  slope  of  the  4  mm 
slit  line  less  than  that  of  the  2  mm  slit  line.  If  a  drastic  cut  in 
secondary  Mixing  length  could  make  the  slope  of  the  4  mm  slit  line  in 
Fig.  32  less  than  the  slope  of  the  2  mm  slit  line,  a  set  of  parameters 
that  accurately  predict  both  2  and  4  mm  slit  total  power  data  could  be 
determined.  The  results  of  the  4  mm  slit  calculation  with  parameter 
Set  D,  shown  in  Table  13  and  Fig.  32,  indicated  no  departure  from  the 
behavior  of  the  previous  4  mm  slit  calculations. 

The  results  of  the  unstable  resonator  calculations  for  the  2  mm  and 
4  mm  slit  widths  indicated  that  one  set  of  baseline  parameters  that 
would  allow  the  model  to  predict  performance  as  the  size  of  the 
resonator  changed  could  not  be  found.  However,  a  set  of  parameters  was 
obtained  for  both  the  2  mm  and  4  mm  slit  cases  that  allowed  the  unstable 
resonator  model  to  accurately  predict  the  total  power  and  oscillation 
period.  The  amplitude  modulation  was  underpred icted  in  the  4  mm  slit 
case.  For  the  2  mm  slit  width,  the  best  fit  parameter  set  was  Set  C, 
Table  12.  For  the  4  mm  slit  width,  the  best  fit  parameter  set  was 
Set  D,  Table  12. 

2.5.1  CL  I  CALCULATIONS 

A  previous  study  demonstrated  that  an  appropriate  set  of  stable 
resonator  baseline  parameters  is  able  to  accurately  predict  laser 
performance  as  the  scale  of  the  laser  device  changes^.  Thus  unstable 
resonator  calculations  were  performed  to  determine  if  parameter  sets  C 
and  D  could  predict  scale  effects.  The  calculations  were  conducted  for 
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the  Hellos  CL  I  laser  for  the  low  pressure.  Run  34  flow  rates.  The 
Hellos  CL  I  laser  Is  Identical  to  the  Helios  CL  II  laser  except  that  the 
flow  rates  and  gain  length  are  reduced  by  a  factor  of  two.  A  2  mm  slit 
calculation  was  performed  with  Set  C  and  a  4  mm  slit  calculation  was 
performed  with  Set  0.  The  results  of  these  calculations  are  shown  in 
Table  14.  Neither  the  2  mm  nor  the  4  mm  slit  calculation  predicted  the 
data.  From  this  It  was  evident  that  parameter  Sets  C  and  D  cannot 
predict  scale  effects. 

2.6  SUMMARY  OF  UNSTABLE  RESONATOR  MODELING  RESULTS 

The  results  of  the  unstable  resonator  baselining  process  indicated 
that  one  set  of  baseline  parameters  that  allow  the  model  to  track 

unstable  resonator  data  as  the  size  of  the  resonator  varied  could  not  be 
found.  It  was,  however,  possible  to  obtain  different  sets  of  baseline 
parameters  for  Individual  resonator  sizes.  The  best  fit  parameter  set 
for  the  2  mm  silt  case  accurately  predicted  the  total  power,  but  also 
predicted  a  20  nsec  small  amplitude  modulation  which  was  not  seen 

experimentally.  The  best  fit  parameter  set  for  the  4  mm  slit  case 

accurately  predicted  the  total  power  and  the  period  of  the  time- 
dependent  oscillations,  but  underpredicted  the  amplitude  modulation. 
Neither  parameter  set  predicted  scale  effects.  The  best  fit  parameter 
set  for  the  5  mm  slit  was  the  FPSR  parameters.  They  predicted  the 

period  and  amplitude  modulation  of  the  oscillations  but  overpredicted 
the  total  power. 

The  Inability  to  obtain  one  set  of  baseline  parameters  for  the 
unstable  resonator  model  demonstrates  that  the  model  is  Inadequate.  The 
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Baseline 

Parameters 


Slit  Total 
Width  Power 
(mm)  (Watts) 


Oscillation 

Period 

(nsec) 


Oscillation 
Aity  1  i  tude 
[%  total  power 


EMM  Data13 

2 

1.97 

— 

— 

EMM  Data13 

4 

0.10 

47 

50 

Set  C 

2 

6.23 

27 

18 

Set  D 

4 

2.48 

27 

36 

Table  14 


Summary  of  2  mm  and  4  mm  slit  MN0R03UR  baselining 
calculations  and  data  for  the  Helios  CL  I  for  the  low 
pressure  Run  34  flow  rates.  All  calculations  used 
tw  =  0.97  and  a  50%  geometric  outcoupled,  confocal, 
unstable  strip  resonator. 


limitations  of  the  model  arise  mainly  from  the  treatment  of  the  fluid 
dynamic  mixing.  Sonic  transverse  Injection  into  a  subsonic  multi- 
component  reacting  shear  flow  results  in  a  three-dimensional  flow 
field.  The  Idealization  of  this  as  a  two-dimensional  array  of  parallel, 
alternating  fuel  and  oxidizer  slit  nozzles  with  subsequent  linear  mixing 
Is  an  oversimplif Ication  which  limits  the  model. 

The  limitations  Imposed  by  the  mixing  model  did  not  prevent  the 
Fabry-Perot  and  stable  resonator  models  from  being  baselined  to  the 
data.  This  suggested  two  possible  explanat ions.  Since  the  actual  fluid 
dynamics  is  independent  of  the  resonator  and  the  unstable  resonator 
model  Is  a  more  exact  treatment  of  the  optics,  the  wave  optics  model  of 
the  unstable  resonator  may  have  brought  out  inadequacies  in  the  mixing 
model  that  were  not  apparent  with  geometric  optics  models.  Another 
possible  explanation  Is  that  the  way  in  which  the  optical  fields  sample 
the  media  in  an  unstable  resonator  has  brought  out  the  inadequa<  'ps  m 
the  mixing  model.  Fabry-Perot  and  stable  resonators  have  no  upstream/ 
downstream  coupling  and  full  upstream/downstream  coupling  of  the  med  a 
and  optical  fields,  respectively.  The  unstable  resonator  coup'es 
upstream  of  the  optical  axis  and  downstream  of  the  optua'  ax's,  but  not 
•cross  the  optical  axis.  Both  possible  explanations  demonstrate  tha* 
unstable  resonator  data  as  a  funit'on  of  resonator  s’/p  and  .  ve  ■;*  the 
laser  provide  a  severe  test  of  a  roup  led,  «a<**  upi :  s,  .  npm  ,  c  _ 


fluid  dynamic  model  of  a  <  hpmx  a 


1  asp' . 


III.  UNSTABLE  RESONATOR  TIME -OE PENCE N!  OSC 1 1 l AT  IONS 


3.1  PROPOSED  MECHANISM  FOR  THE  TIME-DEPENOENT  OSCULATIONS  IN  AN 

UNSTABLE  RESONATOR 

An  objective  of  this  study  was  to  investigate  the  me'^ansm 
responsible  for  the  t ime -dependent  oscillations  that  occur  on  '-res 
whose  saturated  gain  does  not  fill  the  unstable  resonator.  These 
oscillations  have  been  the  subject  of  several  earlier  stud ies ^  • ' '  •  . 

In  these  studies,  the  oscillations  were  found  to  occur  in  both  the  CL  I! 
and  CL  I  lasers  at  low  cavity  pressures  for  resonator  Fresnel  numbers 
above  a  demarcation  Fresnel  number  of  1.5  in  the  Cl  II  and  3.?  in  the 
CL  I.  The  period  of  the  oscillations  was  40  nsec  in  both  the  CL  II  and 
Cl  1  and  was  Independent  of  the  flow  rates.  The  period  of  the  oscilla¬ 
tion  was  determined  by  the  resonator  magnification,  and  the  amplitude 
modulation  was  determined  by  the  fraction  of  the  resonator  filled  by  the 
saturated  gain  of  the  lasing  line.  In  addition  to  the  40  nsec  osc ilia 
tion,  a  7  nsec  osc'liat'on  was  super  imposed  on  it.  The  7  nsec  osc  ilia 
tion  corresponds  to  a  mode  beat  of  the  laser. 

The  data  showed  that  the  occurram.e  of  the  t ’me -dependent  oscilla¬ 
tions  was  Fresne'  number  dependent  and  that  the  period  of  the  oscilla¬ 
tion  was  determined  by  the  resonator  magnification.  Accordingly,  a 
ser'es  of  unstab ’p  resonator  MNOP'  HJP  a'ru'at  vs  with  4  mm  an)  5  mm 
S I  '  t  resonators  were  perfumed  to  ’nvestigatw  the  effeit  of  varying  thp 
resonator  Fresnel  number  and  the  resonator  magnification.  The  resonator 
magnification  was  varied  fhr.jugh  the  radii  of  curvature  of  the  mirrors. 
The  Fresne'  number  was  varied  by  changing  thp  s'/e  of  the  large  mirror. 
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Oc 

4xL 


(3.1-1) 


where  0  Is  the  diameter  of  the  large  mirror,  x  is  the  optical 
wavelength,  and  L  Is  the  mirror  separation.  The  calculations  were 
performed  for  the  CL  II  laser  for  the  low  pressure  Run  34  flow  rates. 
The  baseline  parameters  used  were  the  best  fit  parameter  sets  for  the  4 
mm  and  5  mm  slit  resonators,  Sets  0  and  FPSR,  respectively.  Selected 
Information  from  the  4  mm  and  5  mm  slit  calculations  Is  presented  in 
Table  15. 

The  central  Fresnel  zone  in  an  unstable  resonator  is  the  region 
about  the  optical  axis  whose  radius  Is  such  that  the  Fresnel  number  is 
unity.  In  the  central  Fresnel  zone,  the  media  is  strongly  coupled  to 
the  optical  fields  by  diffraction.  The  size  of  the  central  Fresnel  zone 
Is  determined  by  setting  the  Fresnel  number  to  unity. 


(3.1-2) 


where  0  Is  the  diamter  of  the  central  Fresnel  zone,  L  Is  the  mirror 
separation,  and  x  Is  the  optical  wavelength.  The  diameter  of  the 
central  Fresnel  zone  Is  then 


°CF  =  2/xL 


(3.1-3) 


Thus,  the  extent  of  the  central  Fresnel  zone  is  determined  by  the 
optical  wavelength  and  the  mirror  separation. 
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The  number  of  passes  required  for  a  wave  to  exit  the  resonator 
after  leaving  the  central  Fresnel  zone  is  given  through  a  geometric 


optics  approximation  to  the  unstable  resonator  as 


14 


in 


"p  = 


M(DS) 


12/xl  J 

tn(M) 


(3.1-4) 


where  D$  is  the  diameter  of  the  small  mirror,  M  is  the  resonator 


magnification,  x  is  the  optical  wavelength  and  L  Is  the  mirror 
separation.  Equation  (3.1-4)  shows  that  the  number  of  passes  for  a  wave 
to  exit  the  resonator  is  determined  by  the  resonator  magnification  for  a 
fixed  resonator  size.  It  should  be  noted  that  one  round  trip  through 
the  resonator  corresponds  to  two  passes.  The  number  of  passes  requried 
for  a  wave  to  exit  the  resonator  after  leaving  the  central  Fresnel  zone 
was  calculated  for  the  cases  shown  in  Table  15. 

The  calculations  presented  in  Table  15  showed  a  correlation  between 
the  period  of  the  time-dependent  oscillations  and  the  number  of  passes. 


np,  for  a  wave  to  exit  the  resonator  after  leaving  the  central  Fresnel 


zone.  As  np  Increased,  the  period  of  the  time-dependent  oscillations 
Increased.  The  trend  was  evident  in  both  the  4  mm  and  5  mm  slit 
calculations.  his  correlation  between  the  number  of  passes  for  a  wave 
to  exit  the  re-onator  after  leaving  the  central  Fresnel  zone  and  the 
period  of  the  time-dependent  oscillations  led  to  the  following 
hypothesis  concerning  the  mechanism  responsible  for  the  time-dependent 
oscillations. 


& 


$ 
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The  time-dependent  oscillations  of  the  laser  power  are  the  result 
of  a  competition  between  chemical  pumping  and  radiative  deactivation  due 
to  lasing.  If,  in  a  given  time,  radiative  deactivation  removes  more  HF 
molecules  from  an  upper  laser  level  than  the  chemistry  can  repopulate, 
the  gain  on  that  line  is  reduced.  The  reduced  gain  results  in  reduced 
intensities.  Since  the  stimulated  radiative  processes  are  proportional 
to  the  intensity,  reduced  intensities  reduce  the  radiative  deactivation 
rate.  The  chemical  processes,  which  are  independent  of  the  intensity, 
then  catch  up  and  reestablish  the  population  inversion.  Once  the 
population  inversion  is  reestablished,  the  entire  cycle  repeats 
itself. 

The  data  and  the  MN0R03UR  calculations  in  Table  15  indicated  that 
strong  coupling,  either  diffractive  or  geometric,  prevented  the 
oscillations  by  locking  the  media  and  the  optical  fields  together. 
Strong  diffractive  coupling  occurs  when  the  Fresnel  number  is  small. 
Strong  geometric  coupling  occurs  when  the  photons  require  several  passes 
to  exit  the  resonator.  Thus,  for  the  oscillations  to  occur,  the  rate  of 
chemical  pumping  must  be  less  than  the  rate  of  radiative  deactivation, 
and  the  media  must  not  be  strongly  coupled  to  the  optical  fields  either 
dlff racti vely  or  geometrically. 

3.2  NUMERICAL  INVESTIGATIONS  OF  THE  PROPOSED  MECHANISM 

To  Investigate  the  validity  of  the  proposed  mechanism  for  the  time- 
dependent  oscillations,  several  sets  of  calculations  were  performed.  A 
Fabry-Perot  calculation  was  conducted  with  the  rotational  ronequi 1 ibrium 
MN0R03  model  to  compare  the  relative  magnitudes  of  the  pumping  terms  and 
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the  radiative  deactivation  terms  in  the  species  equation  for  v  =  2. 
These  terms  represent  the  rate  of  HF  production  by  the  pumping  reactions 
and  the  rate  of  HF  depletion  due  to  lasing,  respectively.  The  calcula¬ 
tion  was  carried  out  for  the  low  pressure  Run  34  flow  rates  with  the 
FPSR  baseline  parameters.  The  pumping  and  radiative  deactivation  terms 
from  the  species  equations  for  the  upper  laser  levels  for  two  lasing 
lines  as  calculated  by  MN0R03  are  presented  in  Table  16.  These  results 
indicated  that  the  rate  of  radiative  deactivation  was  two  to  four  times 
faster  than  the  chemical  pumping  rate. 

To  study  the  effect  on  the  time-dependent  oscillations  of  altering 
the  relative  rate  of  chemical  pumping  versus  radiative  deactivation, 
5  mm  slit  MN0R03UR  calculations  using  the  FPSR  parameters  were  performed 
with  the  pumping  reaction  rate  constants  increased.  When  the  pumping 
reaction  rate  constants  were  multiplied  by  a  factor  of  two,  making  them 
gas  kinetic,  the  time-dependent  oscillations  were  slightly  damped.  When 
the  pumping  reaction  rate  constants  were  multipled  by  a  factor  of  ten, 
the  time-dependent  oscillations  were  damped  as  shown  in  Figs.  33  and  34. 

IP 

A  previous  study  had  shown  that  the  time-dependent  oscillations 
occurred  only  when  there  was  rotational  nonequilibrium  present  in  the 
system.  Since  the  concentration  of  HF  Increased  by  a  factor  cf  ten  «he’ 
the  pumping  rate  was  Increased  by  a  factor  of  ten  and  the  -ota^ora" 
relaxation  rate  increases  as  the  square  of  the  HF  concent-at  0"  ,3 
factor  of  one-hundred  in  cis  case),  an  MN0Ru3uR  caicu'at  5on 
performed  to  Investigate  whether  the  damping  that  occurred  t*e 

pumping  rate  constants  were  multiplied  by  ten  was  due  te  an  increased 
rotational  relaxation  rate.  The  pumping  rate  constants  were  set  at 
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their  experimental  values  and  the  rotational  relaxation  rate  was 
increased  by  a  factor  of  one  hundred.  Figures  35  and  36  show  that,  with 
the  rotational  relaxation  rate  increased  by  factor  of  one  hundred,  the 
amplitude  modulation  of  the  time-dependent  oscillations  decreased  and 
the  period  was  unchanged  compared  to  the  a  priori  calculation  (Figs.  3 
and  4).  The  damping,  however,  was  not  as  great  as  when  the  pumping  rate 
constants  were  increased  by  a  factor  of  ten.  Figs.  33  and  34.  This 
indicated  that  the  increased  rotational  relaxation  rate  was  not 
responsible  for  the  heavy  damping  seen  in  the  calculation  with  the 
pumping  rate  constants  multiplied  by  a  factor  of  ten. 

The  preceding  two  MN0R03UR  calculations  showed  that  increasing  the 
chemical  pumping  rate  damped  the  time-dependent  oscillations,  and  this 
damping  was  due  to  the  increased  rate  of  production  of  HF  molecules  by 
the  pumping  reactions,  not  the  increased  rotational  relaxation  rate. 
These  results  supported  the  hypothesis  that  the  mechanism  responsible 
for  the  time-dependent  oscillations  is  a  competition  between  chemical 
pumping  and  radiative  deactivation  of  the  upper  laser  levels  of  HF. 

According  to  the  proposed  mechanism,  if  the  media  is  strongly 
coupled  to  the  optical  fields,  the  oscillations  are  damped  out  since  the 
media  and  the  optical  fields  are  locked  together.  Previous  MN0P03UP 
calculations,  Table  15,  Illustrated  the  effect  on  the  time  -dependent 
oscillations  of  geometric  coupling  of  the  med 1  a  an.:  op*'  a'  f  is.  fcr 
a  fixed  resonator  Fresnel  number,  as  np  ’ni.reaspd,  per-,,,<  jf  the 
time -dependent  oscillations  increased. 

Oiffractive  coupling  of  the  med’a  anil  *  nP  r*.;  ,n  *'p\j  «as 
investigated  with  mnoP'jIUP  ra  «•  j  'at  ions  ’.ra’  «a'  **  •  •  np  rp.  .rv  >■ 


Outcoupled  Power,  watts 


HELIOS  CL  H  MN0R03UR 

Run  34  low  pressure 

50%  geometric  outcoupled 

Confocal  unstable  resonator 

NF=  8  93  0 A  =  0  5cm 

T,«  =  T|ws  10 

FPSR  parameters 

Rotational  Relaxation  Rote  x  1C 


Iteration  Number 


Figure  35  MNOR03UR  multiline  outcoupled  power  calculated  -n  the 
optics8,  PQ  ^  ,  vs.  iteration  number  for  rota*'una 
nonequ  i  1  i  brTum  for  a  bO%  geometric  out<oup'ed,  <oMu<<i 
unstable  strip  resonator  with  a  large  m'rror  o  i  arm- '  *- * 
1.0  cm  for  the  Helios  Q  !I  laser  for  Pur  34  at  t .  4 


or.  .  . 
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Fresnel  number  by  changing  the  mirror  spacing.  The  calculations  were 
performed  for  the  CL  II  laser  with  a  4  mm  slit  resonator  for  the  low 
pressure.  Run  34  flow  rates.  Baseline  parameter  Set  0,  Table  12,  was 
used.  Table  17  shows  three  MN0R03UR  calculations.  In  calculation  A, 
the  nominal  calculation,  the  mirror  spacing  was  100  cm  and  the  resonator 
Fresnel  number  was  5.71.  In  calculation  B,  the  mirror  spacing  was  200 
cm,  double  the  nominal  value,  and  the  resonator  Fresnel  number  was 

2.85.  In  calculation  C,  the  mirror  spacing  was  50  cm,  half  the  nominal 
value,  and  the  resonator  Fresnel  number  was  11.42.  Calculation  A,  with 
a  Fresnei  number  of  5.71,  exhibited  40  nsec  time-dependent  oscillations 
as  seen  in  Figs.  37  and  38.  Calculation  B,  with  a  Fresnel  number  of 

2.85,  had  a  central  Fresnel  zone  that  was  larger  than  the  small 
mirror.  Thus,  a  large  portion  of  the  active  media  was  in  a  region  of 
high  diffractive  coupling.  Calculation  B  showed  no  time-dependent 
oscillations.  The  results  of  calculation  C,  with  a  Fresnel  number  of 
11.42  are  presented  in  Figs.  39  and  40.  A  comparison  of  the  results  of 
calculations  A  and  C  shows  that  as  the  Fresnel  number  was  increased,  the 
amplitude  of  the  time-dependent  oscillations  increased.  This  was  in 
accordance  with  the  proposed  mechanism,  which  indicates  that  as  the 
degree  of  diffractive  coupling  between  the  media  and  the  optical  fields 
decreases,  the  oscillations  become  stronger.  In  this  case,  the 
diffractive  coupling  decreased  (larger  Fresnel  number)  while  the 
geometric  coupling  stayed  the  same  (np  approximately  the  same  as  for 
case  A).  In  addition,  the  results  of  calculation  C  showed  that  the 
period  of  the  oscillations  decreased.  This  was  a  result  of  the 
reduction  In  round  trip  transit  time  due  to  the  decreased  mirror 
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gure  37  MN0R03UR  multiline  outcoupled  power  calculated  in  the 
optics8,  PQ(jt,  vs.  iteration  number  for  rotational 
nonequilibrium  for  a  50%  geometric  outcoupled,  confocal 
unstable  strip  resonator  with  a  large  mirror  diameter  of 
0.8  cm  for  the  Helios  CL  II  laser  for  Run  34  at  5.4  Torr 


tcoupled  Power,  watts 


Outcoupled  Power,  watts 
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Figure  39  MN0R03UR  outcoupled  power®,  PQU^,  for  individual  lines  vs. 

iteration  number  for  a  50%  geometric  outcoupled,  confocal, 
unstable  strip  resonator  with  a  large  mirror  diameter  of 
0.8  cm  for  the  Helios  CL  II  laser  for  Run  34  at  5.4  Torr. 


Outcoupled  Power,  watts 


The  data  showed  that  the  oscillations  occurred  for  Fresnel  numbers 
above  1.5  in  the  CL  II  and  3.2  in  the  CL  I.  This  difference  in  the 
demarcation  Fresnel  numbers  is  due  to  the  gain  length  in  the  CL  I  being 
half  that  in  the  CL  II.  Because  the  gain  length  in  the  CL  I  is  shorter 
and  the  diffractive  losses  are  higher,  for  the  same  resonator  Fresnel 
number  (Eq.  (3.1-1)),  the  CL  I  operates  at  an  effective  Fresnel  number 
that  is  lower.  Thus  for  the  oscillations  to  appear,  the  resonator 
Fresnel  number  must  be  higher  in  the  CL  I  than  in  the  CL  II. 

3.3  RELATIONSHIP  OF  THE  PROPOSED  MECHANISM  TO  OTHER  MODELS 

Other  models  of  the  temporal  stability  of  unstable  resonators  have 
been  proposed  by  Dreizin  and  Dykhne^  and  Mirels^®.  The  present 
hypothesis  considers  the  time-dependent  oscillations  as  the  result  of  a 
competition  between  chemical  pumping  and  radiative  deactivation  of  the 
upper  laser  levels  of  HF.  The  models  of  Dreizin  and  Dykhne  and  Mirels 
consider  the  time-dependent  oscillations  as  the  result  of  a  chemical 
kinetic-fluid  dynamic  phenomena,  in  which  fluid  elements  flow  into  the 
upstream  portion  of  the  resontor,  lase,  and  are  subsequently  replaced  by 
fresh  media  in  the  upstream  portion  of  the  resonator.  Accordingly,  the 
oscillation  period  is  on  the  order  of  the  fluid  transit  time  through  the 
resonator.  The  experimental  evidence^’ ^ however,  indicates  that 
time-dependent  oscillations  occur  even  when  the  velocity  of  the  media  is 
so  low  that  the  media  is  essentially  frozen  for  times  on  the  order  of 
the  period  of  the  oscillations.  Therefore,  it  is  clear  that  these  time- 
dependent  oscillations  are  not  the  result  of  chemical  kinetic-fluid 
dynamic  coupling.  The  mechanism  responsible  for  these  time-dependent 


oscillations  Is  a  competition  between  chemical  pumping  and  rad'at  we 
deactivation. 

3.4  SUGGESTED  EXPERIMENTS 


The  proposed  mechanism  for  the  time-dependent  oscillations  in 
unstable  resonators  on  lines  whose  saturated  gain  does  not  fill  the 
resonator  was  supported  by  the  calculations  in  Section  3.2.  To  verify 
the  proposed  mechanism,  new  unstable  resonator  experiments  should  be 
conducted.  Based  on  the  proposed  mechanism,  for  a  fixed  resonator 
Fresnel  number,  the  period  of  the  time-dependent  oscillations  should 
increase  as  the  geometric  coupling  increases,  since  the  media  is  more 
coupled  to  the  optical  fields.  Thus  an  experiment  with  a  4  mm  slit 
resonator  with  a  magnification  of  1.25  (20%  outcoupled)  should  be 
performed  to  determine  whether  the  period  increases  as  predicted.  The 
proposed  mechanism  also  predicts  that  the  oscillation  amplitude 
modulation  should  Increase  as  Fresnel  number  increases  since  the  media 
is  less  coupled  to  the  optical  fields.  If  the  geometric  coupling  is 
held  constant  as  the  Fresnel  number  is  increased,  the  period  should 
decrease  as  shown  by  the  MN0R03UR  calculations  in  Tables  15  and  17. 
Thus  experiments  with  mirror  spacings  of  200  cm  and  50  cm  should  be 
performed  for  a  50%  geometric  outcoupled  4mm  slit  resonator.  There 
should  be  no  oscillations  observed  with  the  200  cm  mirror  spacing  since 
the  diffractively  coupled  central  Fresnel  zone  will  be  larger  than  the 
small  mirror  and  contain  most  of  the  media.  The  oscillations  observed 
with  the  50  cm  mirror  spacing  should  have  a  larger  amplitude  and  a 
shorter  period  than  those  for  a  ’00  cm  mirror  spacing  since  the  media 
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»'for  spac’ng  snou’d  a  ’  sc  change  t*>e  pe^od  of  the  node  Beat,  for  a 
200  cm  mirror  spacing,  the  mode  Beat  shou’d  have  a  period  of  13  tc 
14  nsec.  For  a  50  cm  mirror  spacing,  the  mode  Beat  should  not  appear 
since  the  longitudinal  mode  frequency  of  approximately  300  MHz  and  the 
line  width  of  approximately  400  MHz  indicate  that  only  one  longitudinal 
mode  would  lase. 

3.5  SUMMARY 

A  mechanism  to  explain  the  time-dependent  oscillations  that  occur 
on  lines  whose  saturated  gain  does  not  fill  the  unstable  resonator  was 
proposed.  It  was  suggested  that  the  time-dependent  oscillations  are  the 
result  of  a  competition  between  chemical  pumping  and  radiative 
depopulation  of  the  upper  laser  levels  of  HF.  The  oscillations  occur 
only  if  the  media  is  not  strongly  coupled  to  the  optical  fields 
diffractively  or  geometrically. 

Computer  calculations  with  the  full,  nonlinear,  coupled,  chemical 
kinetic,  fluid  dynamic,  wave  optics  unstable  resonator  model  supported 
the  proposed  mechanism.  The  calculations  showed  that  the  oscillations 
were  damped  out  by  either  increasing  the  rate  of  the  chemical  pumping 
reactions  or  by  increasing  the  degree  of  coupling  between  the  media  and 
the  optical  fields. 

Three  experiments  with  a  4  mm  slit  resonator  were  suggested  to  test 
the  validity  of  the  proposed  mechanism.  For  a  fixed  resonator  Fresnel 
number,  decreasing  the  magnification  of  the  resonator  to  1.25  should 
increase  the  period  of  the  time-dependent  oscillations  since  the  media 
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will  be  more  geometrically  coupled  to  the  optical  fields.  Changing  the 
mirror  spacing  to  50  cm  should  increase  the  amplitude  modulation  and 
decrease  the  period  since  the  media  will  be  less  diffractively  coupled 
to  the  optical  fields.  Changing  the  mirror  spacing  to  200  cm  should 
prevent  the  oscillations  since  the  central  Fresnel  zone  will  be  larger 
than  the  small  mirror  diameter,  thereby  increasing  the  diffractive 
coupling  of  the  media  to  the  optical  fields. 
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IV.  CONCLUDING  REMARKS 

The  objective  of  this  study  was  to  simulate  cw  HF  chemical  laser 
unstable  resonator  performance.  To  accomplish  this  the  unstable 
resonator  data  were  modeled  with  a  fluid  dynamic,  chemical  kinetic, 
rotational  nonequilibrium,  wave  optics  model,  MN0R03UR.  The 
calculations  showed  that  one  set  of  baseline  parameters  that  allow  the 
model  to  predict  the  data  as  the  size  of  the  resonator  and  size  of  the 
laser  varied  could  not  be  found.  This  inability  to  obtain  one  set  of 
baseline  parameters  revealed  limitations  in  the  unstable  resonator 
model.  The  limitations  are  due  primarily  to  the  treatment  of  the  fluid 
dynamic  mixing.  Since  these  limitations  were  not  apparent  in  Fabry- 
Perot  and  stable  resonator  models  with  the  same  mixing  model,  the 
unstable  resonator  data  as  a  function  of  resonator  size  and  size  of  the 
laser  was  shown  to  be  a  severe  test  of  a  coupled,  wave  optics,  chemical 
kinetic,  fluid  dynamic  model  of  a  chemical  laser. 

To  improve  the  MN0R03UR  model,  the  treatment  of  the  fluid  dynamic 
mixing  should  be  Improved.  Sonic  transverse  injection  into  a  subsonic 
multicomponent  reacting  shear  flow  results  in  a  three-dimensional  flow 
field.  The  Idealization  of  this  as  a  two-dimensional  array  of  parallel, 
alternating  fuel  and  oxidizer  slit  nozzles  with  subsequent  linear  mixing 
Is  an  oversimplification.  A  full,  three-dimensional  mixing  model  is 
needed. 

A  mechanism  that  explains  the  time-dependent  oscillations  that 
occur  on  lines  whose  saturated  gain  does  not  fill  the  unstable  resonator 


was  proposed.  According  to  the  proposed  mechanism,  the  time-dependent 
oscillations  are  the  result  of  a  competition  between  chemical  pumping 
and  radiative  deactivation  of  the  upper  laser  levels  of  HF.  The 
oscillations  occur  only  If  the  media  Is  not  strongly  coupled  to  the 
optical  fields  dlffractlvely  or  geometrically.  Computer  calculations 
supported  the  proposed  mechanism.  To  verify  the  proposed  mechanism,  new 
unstable  resonator  experiments  should  be  conducted.  Based  on  the 
proposed  mechanism,  for  a  fixed  resonator  Fresnel  number,  the  period  of 
the  time-dependent  oscillations  should  Increase  as  the  magnification  of 
the  resonator  decreases.  Thus  an  experiment  with  a  4  mm  slit  resonator 
with  a  magnification  of  1.25  (20%  outcoupled)  should  be  conducted.  For 
this  resonator,  the  period  of  the  oscillations  should  be  greater  than 
for  a  4  mm  slit  resonator  with  a  magnification  of  2.  The  proposed 
mechanism  also  predicts  that  the  oscillation  amplitude  modulation  should 
increase  and  the  period  decrease  as  the  Fresnel  number  of  the  resonator 
Increases  with  the  geometric  coupling  held  constant.  Thus  experiments 
with  mirror  spaclngs  of  200  cm  and  50  cm  should  be  performed  for  a  50% 
geometric  outcoupled  4  mm  slit  resonator.  There  should  be  no  oscilla¬ 
tions  with  the  200  cm  mirror  spacing,  except  for  the  mode  beat  of  the 
laser  which  should  have  a  period  of  13  to  14  nsec.  With  the  50  cm 
mirror  spacing,  there  should  be  oscillations  that  have  a  larger 
amplitude  and  a  shorter  period  than  the  oscillations  seen  with  a  mirror 
separation  of  100  cm.  A  mode  beat  with  the  50  cm  mirror  spacing  should 
not  be  observed. 
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